Instability and decomposition on the surface of strained alloy films 
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A continuum dynamical model is developed to determine the morphological and compositional 
instabilities on the free surface of heteroepitaxial alloy films in the absence of growth. We use linear 
stability analysis to study the early nonequilibrium processes of surface evolution, and calculate 
the stability conditions and diagrams for different cases of material parameters. There are two key 
considerations in our treatment: the coupling between top free surface of the film and the bulk phase 
underneath, and the dependence of both Young's and shear elastic moduli on local composition. 
The combination and interplay of different elastic effects caused by lattice misfit between film and 
C^) ' substrate (misfit strain), composition dependence of film lattice parameter (compositional strain) 

and of film elastic constants lead to complicated and rich stability results, in particular the joint 
stability or instability for morphological and compositional profiles, the asymmetry between tensile 
and compressive layers, as well as the possible stabilization and suppression of surface decomposition 
even below the effective critical temperature. We also compare our results with the observations of 



. some postdeposition annealing experiments. 

I. INTRODUCTION 
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Instability is one of the most important phenomena in strained heteroepitaxial growth, where a thin film of different 
material is epitaxially deposited on a substrate and the mismatch of lattice constants between film and substrate leads 
to strain in the grown film. Besides the formation of misfit dislocations and other defects, the strain can be relieved 
I | through the process of morphological instability, during which the growth mode of the film is changed from layer- 
' O i by-layer to three-dimensional island but the film still remains coherent with its substrate. This dislocation-free 
morphologicaLinstability occurs in a|-.wide range of heteroepitaxial systems, and has been well investigated both 
theoreticallylaffiicl and experimentally0i3 for single-component materials . 

In recettUvftarSyimore attention has been paid to the mu4t| VCQmponent strained layers, in particular the binary (e.g., 
Sii-KGe^Lrffl'td'EJ) and pseudo-binary (e.g., Ini_ K Ga x As^'^oE3ll30) alloys. Since the materials are composed of 
different kinds of atoms that may not be fully miscible, the compositional inhomogcncitics arc expected to develop at 
certain growth conditions. This alloy segregation instability couples with the morphological instability msuMng in the 
■ simultaneous modulations of the surface profile and alloy composition during the epitaxial growth. BE3ESli3 Another 
important phenomenon of the heteroepitaxial alloy films is the asymmetry of the stability for compressive andptensile 
, layers. For Sii-^Ge^, films grown under compression were observed to be less stable than those under tension,0 while 
for Ini_ x Ga x As the situation is more-cpi»Blicated: Whether compressive or tensile films are more stable seems to 
, depend on the material's growth rate.EjIi-ftj Compared with single=£jomponent films, there is an extra effect on alloy 
strain due to composition dependence of the film lattice constant S3i3 The combination f^t-has compositional strain 
and the film-substrate misfit strain highly influences the stability behaviors of alloy films£3'Eil 

Some theoretical investigations have been carried out in order to understand the morphological instability of strained 
miilti-pnmppnprit films and most of them focused on the properties and behaviop-xdj growing layers with deposition 
flux.t3EJEjOoE3cj However, the studies of the static film without depositionoO are much fewer, although the 
stability of epitaxial films in the absence of growth is also important from both the theoretical and experimental 
points of view. Using the treatment of thermodynamics, GlasEO has demonstrated that any strained alloy with a 
O ■ static free surface is unstable due to the interplay of morphological and compositional instabilities. Furthermore, 
as also pointed out by Glas, this result is valid provided all the evolution mechanisms are possible, and when the 
ijji ' physically relevant mechanism is considered, the unstable state derived thermodynamically may not be, kinetically 
accessible. This has been verified in the nonequilibrium and dynamical analysis of Leonard and Desai,EJ where the 
stability properties of both the non-growing (static) and growing strained alloy films were determined. Stabilization 
in non-growing film was found to be possible for some parameter values, showing the significance of kinetic evolution 
process. 

In this paper we apply a nonequilibrium, continuum model in order to study the morphological instability and 
surface decomposition of binary or pseudo-binary strained non-growing alloy films, i.e., films without deposition. 
The system is assumed to be clastically isotropic and follow the conserved dynamics of phase separation and surface 
diffusion. Compared with the model of Leonard and Desaicfl and other previous work, here we have two crucial 
considerations which lead to significantly different stability results. First, due to the fact that the surface phase and 
bulk phase in a strained film are intimately coupled with each other, we take into account the total free energy of 
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the system, not just the surface state and energy as used beforeS Second, the film elastic moduli (both the Young's 
modulus E and shear modulus fi) are considered to be composition-dependent, which deeply affects the stability 
properties. Presenting the linear analysis on early evolution of surface profile and composition, we derive the stability 
conditions for non-growing film as well as stability diagrams for various material parameters appropriate to realistic 
systems. 

In Sec. II, wc describe the continuum model with conserved dynamics and present the mechanical equilibrium 



equation for the elasticity of this film system. The linear stability analysis is carried out in Sec. Ill, resulting in the 
characteristic equation for perturbation growth rate. The results related to the stability conditions and diagrams are 
shown in Sec. IV, and the discussion and conclusions are presented in Sec. 



II. MODEL AND ELASTICITY FORMULATION 

We consider a strained alloy system composed of a semi-infinite substrate occupying the region z < and a h\-x^>x 
binary or pseudo-binary alloy film in the region < z < h(x,y), where h(x,y) represents the surface height variable. 
The film-substrate interface located at z — is assumed to be planar and remain coherent without generating any 
misfit dislocation. In this hcteroepitaxial system, the misfit strain caused by the difference between the lattice constant 
of the film af and that of the substrate a s is characterized by e = (a/ — a s )/a s . Thus, e > or < implies a strained 
film under compression or tension, respectively. To describe the composition profile of the film, we use a continuous 
variable </>(r, t) which is defined only for z > and is proportional to the local difference in the concentrations of two 
constituents. Corresponding to an alloy composition X, its average value <j) is equal to 2X — 1. Here we focus on 
the symmetric mixture, i.e., X — 1/2 alloy, for which the spinodal decomposition theory can be well applied, and 
then we have cf> = in what follows. Due to the atomic size difference between A and B,|Seecies, the film lattice 
parameter af is composition-dependent and the solute expansion coefficient rj = d In af /d(j&^3 is defined to measure 
the compositional strain. 

Assume that the thin film has been grown under ultra-high vacuum condition, e.g., in a molecular-beam epitaxy 
(MBE) system, and then the evaporation and re-condensation on film surface are negligible. Furthermore, we also 
neglect the interdiffusion between film and substrate as well as the diffusion and compositional relaxation in the bulk 
film, since the bulk atomic mobility is much smaller than the mobility at the surface in typical epitaxial growth. 
Thus, the dynamics of morphological and compositional evolution is dominated by the surface diffusion and surface 
decomposition processes, and should be conserved. 

For the evolution of surface profile, the surface diffusion mechanism leads to 

while to measure the time-dependence of concentration field at the surface 4>{x, y, h(x, y),t)— <f> s {x, y, t), we apply the 
conserved dynamics: 

Note that in Eq. (|^) there are two ways to study the surface composition fluctuations. The first oneS is evaluating 
the free energy T at the surface and then calculating the functional differentiation with respect to surface composition 
field 4> s , that is, only considering the surface .state. Here we take into account the intimate coupling between surface 
state and bulk state and use the other way:E3 Apply the total free energy T of the whole system to calculate the 
composition dynamics and then evaluate it at the surface, as has been done in the previous study of surface critical 
phenomena for spin fluctuations .Ej 

In Eqs. (Q) and (Q), V 2 is the surfaoe-Laplacian, g = 1+ |V/i| 2 represents the determinant of the surface metric, and 
the kinetic coefficients are denoted asc2l Th = D s N s /kBTN% and = T^S -1 with S the effective diffusion thickness 
of surface layer. Here D s is the surface diffusivity, fc^ is the Boltzmann constant, N s and N v are the number densities 
of atoms per unit surface area and per unit volume, respectively, and T is the temperature. The total free energy 
functional T consists of three contributions: 

F = F s + F Gh + F eh (3) 

The first contribution T s is the surface energy, which plays a stabilizing role and can be represented by a drumhead 
model without pinning term: 



TM=1 J d 2 r^g. (4) 
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Here 7 is the surface tension, and for simplicity we assume it to be isotropic and composition independent. The 
second term in r.h.s. of Eq. (||) determines the phase behaviors of binary compounds and is the Ginzburg-Landau 
functional 



F Gh [<f>,h] = / d 3 

J —OO 



f|V*|* 



(5) 



where the parameters^! r' = fcs(T c — T)N V , u is a temperature independent positive constant, and c = kBT c N v a^/2, 
with T c the critical temperature of the binary alloy and do the effective interaction distance. For the bulk alloy without 
elastic strain, when T > T c the equilibrium state is homogeneous with <fi = 0, while for T < T c we have the coexistence 
of two phases <fi = ±^/r' /u. The last term |V^| 2 in Eq. (^) represents the gradient energy that penalizes the sharp 
compositional changes, rand, is important for stability analysis: The lack of it leads to a nonphysical divergence for 
short wavelength modeEJE3 

The last contribution in Eq. (|^) is the elastic free energy functional T e i, and is crucial for this stress-driven system. 
From linear elasticity theory, it can be expressed as 



1 f h 



(6) 



where is the stress tensor and Sijki is the elastic compliance tensor with the form Sijki — 5ikSji(l+u)/ E — bijE^iv j E 
for isotropic systems (subscripts i, j, fc, or I = x, y, z). Generally, the elastic constants (Young's modulus E, shear 
modulus n, and Poisson ratio v) are dependent of the local composition, and here we consider this dependence to first 
order, that is, 



E = Eo(l 

fi = fi (i + n\ 



(7) 



and v = E/2fi — 1. In this paper we take the system as being elastically isotropic, and neglect the difference in the 
average elastic constants (Eq, fiQ, and vq = Eq/2^jlq — 1) between film and substrate. This is appropriate for systems 
with substrate and film having similar elastic constants. 

To determine the elastic energy (^|), we need to get the solution for the displacement vector u which satisfies 
mechanical equilibrium 







(8) 



in the whole film/substrate system. According to Hooke's law for isotropic system, the linear stress tensor is expressed 
as 



2/z 



1 



-(e + r)<j>)5i. 



1 -2v " J J \-2v 

with the presence of misfit strain e and composition strain rj<p, where the linear strain tensor mj is given by 

Uij = (diUj + djU.i)/2. 



(9) 



(10) 



The boundary conditions are needed to solve the above mechanical equilibrium equation. At the free surface of the 
film, i.e., at z = h(x,y), we have 







(11) 



due to the negligible pressure on the film surface. Here rij is the unit vector normal to the surface. Since the film- 
substrate interface at z = remains coherent, we get the continuous conditions for both stress and displacement 
tensors: 



4 = < 3 



and 



u{ = uf 



(12) 



where superscripts / and s refer to the film and substrate, respectively. Finally, the strains far from the interface, 
that is, for z — > —00, are expected to decay to zero: 







and 



0. 



(13) 
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III. LINEAR STABILITY ANALYSIS 



In order to determine the stability properties of this non-growing strained alloy system, we apply the linear analysis 
on evolution equations (|j) and (g) with the use of formulae (pf)-(|l3|). For a general variable £, which could be 
displacement itj, composition field </>, or height variable h, its Fourier expansion yields 



q 



i(q,z,t)e 



i(q^x+q y y) 



(14) 



with small perturbations £ around the basic state £ which corresponds to a planar film with fixed thickness ho = h 
and uniform composition_^_p 0. Note that in Eq. (14) when £ denotes the height variable h, £ is in fact /i(q, i). 



The basic-state solutionoo for the film leads to u 



-J 



if 



if 



-2{1qu, and other stress and strain tensors {u{, 



is unstrained and then 



and 



0, u{ 
J 



uz with u = u( z = e(l + vq)/(1 — vq), 
®It) t° be zero. For the substrate, the basic state 

~4j = d h = (*) 3 = y, z). 

The mechanical equilibrium equation (JsJ) with boundary conditions (pd|)-(p^|) can be solved toJirst order 0(h, (j>) 
using the above expansion (U). Here we use the detailed solutions given by Leonard and DesaO, where a crucial 
step is to introduce a new variable W with 



V 2 W 



(15) 



or equivalently, (d% — q 2 )W — with q 2 = q 2 + ql- After substituting the solutions in the free energy functional (||)-(||) 
and then in the dynamical equations (Q) and (^|), we can obtain the linearized evolution equations for morphological 
and compositional perturbations h and <fi to determine the stability of the system. 

What we are interested in is the behavior of perturbations for the stressed film without deposition. The fluctuations 
of alloy composition mainly occur at the surface z — h(x, y) due to the surface relaxation process, and should attenuate 
along the vertical direction z as the surface/bulk coupling weakens with the increasing distance from the surface and 
the bulk mobility is very small. Thus, the bulk compositional perturbation <pb caused by the free surface disturbance 
is hypothesized to decay as 



-k{Ho— z) 



(16) 



with the corresponding W — 4> s e^ K ^Z^ /(n 2 — q 2 ). This exponential form is for the early evolution regime and 
similar to that used in previous work.otll The parameter k in Eq. ( |l6| ) is equal to 1/6, with b the vertical length scale 
of compositional perturbation caused by free surface. Due to the negligible atomic mobility in the bulk and guided 
by the fact that the vertical morphological perturbation is very small compared with the lateral variation, one can 
assume b -C A, where A ~ l/q is the typical lateral wavelength of surface modulation. Therefore, for the range of q 
that corresponds to typical surface structure, we have 



« » q 



(17) 



in Eq. @. 

Using the solutions of mechanical equilibrium equation and the assumptions ( |l6| ) and ([l7|), we have derived the 
dynamical equations for h and <p s . To first order of the perturbations, they are (in nondimensional form) 



dh*/dr= {e* 2 k 3 - 7 *fc 4 )/i 



l + ^o 



2£ 1 '-(1 + ,qK ,a 
V 2(1 -i*>) 



(18) 



and 



l/dr = 



(i - 2is y v * + (2£* - (i + vo)fiy 



8^-5(1 + ,qK 
+ 2(1 - v 2 ) 6 71 



(19) 



where the "±" sign corresponds to the cases in which the alloy is above ("+") or below ("— ") the effective critical 
temperature T® ff defined by 



T cff = T 



2E V 2 
1-vq ksN v 



(20) 
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which is the same as the spinodal temperature of alloys with coherency stress derived by Cahn.0 

Here we have rescaled the variables and parameters to make the equations nondimcnsional, using a characteristic 
length scale 



-1/2 



which is the typical width of domain interfaces, and time scale 



to = ( r^— - 



with 



as well as the transformations: 



r = r' - 2E rf/{l - vq) = M^f - T)N V , 



(21) 



(22) 



(23) 



V = 









t/r , 






h/l a , 












lo 

—1, 

c 






~2E 


{t 


+ vo 


\r\ 




- v a 


~2E 


(f 


+ "0 


. \r\ 




- Vq 



1/2 



-,1/2 



(24) 



In the early time regime of the perturbation's evolution, the growth rates of morphological and compositional 
perturbations, i.e., Oh and o^, are defined through h* = ho exp^r) and d>* = </>o ex p(°> r )! respectively. In general 
cases (e.g., both e* and rf are nonzero) two dynamical equations ( |18| ) and (|19| ) couple with each other, which leads to 
the joint stability or instability of surface morphology and composition, that is, <j] x = = a. Therefore, from Eqs. 
( |l8| ) and (|l^) we can obtain the characteristic equation for the perturbation growth rate er: 



(a + 7*fc' - e 
k 5 



k 2 k 2 ± 1 



2{l + v ){l-v {) f 



8EI - 5(1 + , . 

2(1 -*g) £7? 

(1 - 2v )e*r)* + (2^1* - (1 + ^oKK 



x 2(1 - iy )eV + (2^i* - (1 + ^oKK =0 



(25) 



where the "+" ("— ") sign corresponds to T > T° s (T < T° ). The real part of a determines the stability properties 
of the system: For Re(cr) > 0, the film is jointly unstable to morphological instability and alloy decomposition 
instability at the surface, while if Rc((t) < is fulfilled for all the wavenumbers k, both the morphological and 
compositional modulations are supp resse d. When Re (a) > and Im(cr) ^ 0, the instability is oscillatory. In the 
following calculations, we first (Sec. IV A and IV B) focus on the real part of a for each of the solutions of Eq. (|2^), 
and use the largest one (with respect to all evolution mod es fc) to determine the regions of instability. Within the 
unstable region, the imaginary part is then computed (Sec. IV C) to determine the regions of oscillatory instability. 

Note that in Eq. (p5|), quantities k, 7*, e* , if depend on temperature via the coefficient r (see Eqs. ( p3f), (pi]), 
and (pi])). The temperature dependence of the system's stability can be made explicit by rewriting Eq. (|25|) in a 
dimension-full manner, which is 



c 3 / 2 



cq 



l-u 
E 



(i - 



r (8£S-5(l + n,K)o7 
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g 5 [(1 - 2^ )«7 + (2££ - (1 + n>K)e 2 



M 9 / 2 (l-^o) 4 

x [2(1 - v )er) + (2E{ - (1 + v )iA)<?] = 0- (26) 

Thus, in this characteristic equation for the perturbation growth rate a (= £Itq), the temperature dependence enters 
only through the quantity r which is linearly dependent on T as given by Eq. (|2jj). Even though the characteristic 
time scale To depends on temperature (Eq. (53)), this dependence does not have any effect on stability boundaries. 

IV. RESULTS 

From Eqs. (|l8|)-(^5|), the stability property of the strained film depends on material parameters e*, ij* , E*, fj,*, 7*, 
and vq. In the following we give the results of film stability for the cases of composition independent and dependent 
clastic constants, as well as different conditions of misfit and compositional strains. 

A. Composition-independent elastic moduli (E* = fj,* = 0) 

When ignoring the dependence of elastic constants on the local composition, we have E = Eq, fi = fi$, and v = vq 
from Eq. (Q), and the derived evolution equations for perturbations h* and (jf s as well as the characteristic equation 
for a are the same as Eqs. @, @ and @ after setting E{ = fj,* = 0. 

When e* = rj* = 0, that is, neither misfit strain nor compositional strain exists in the film, the dynamical equations 
for h* and <j>* decouple, with different perturbation growth rate: 

o h = - 7 *fc 4 , 

_ f _ fc 4 _ fc2) if T > Tf 
a * ~ \ -k A + k\ if T < Tf [27) 

which recovers the results obtained by Leonard and Desai.0 Thus, the surface morphology is always stable and the 
compositional stability is similar to that of bulk alloy: For T > T° ff (= T c , when r\ = 0) the system is stable, while 
for T < T° s spinodal decomposition occurs for long wavelengths (k < 1). For the case of zero misfit but nonzero 



solute expansion coefficient, i.e., e* = and rj* ^ 0, the dispersion relations are the same as (27). Note that the 
compositional perturbation rate er^ obtained here is different from that of the previous work (see Eq. (44) in Ref. 
24). 

When e* 7^ and rf = 0, corresponding to nonzero misfit stress but zero solute stress, our results yield: 

o h = e * 2 k 3 -j*k\ 

-fc 4 - fc 2 , if T > T c 

-k A + fc 2 , if T<T r [Z *> 



where the Asaro-Tiller-Grinfcld instability for morphology^ is recovered, as also shown in the model of Leonard and 
Desai. 

For more general case of e* 7^ and 77* 7^ 0, that is, the strains generated by both the lattice mismatch and 
compositional nonuniformity are nonzero and coupled, we have a quadratic equation for the common perturbation 
growth rate a: 

a 2 + a X (T + a = 0, (29) 
with the roots a = a\ ± \J a\ — 4ao^ /2 and coefficients 

ai = 7 *fc 4 -e* 2 fc 3 + fc 2 (fc 2 ±l), 

a - fc 2 (fc 2 ±l)( 7 *fc 4 -e* 2 fc 3 ) + i-^e*V 2 fc 5 . (30) 

1 v Q 

Usually the Poisson ratio vq lies in the range from 1/4 to 1/3. Consequently, for T < Tf (bottom sign "— " in 
Eq. (p0|)), we have a,\ < and clq > when the wavenumber is very small, i.e., k 1, corresponding to the 
solution Re(cr) > 0. Then the instabilities of surface morphological and compositional profiles are expected to appear 
simultaneously below the effective critical temperature. 
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For T > T? ff , the stability properties are more complicated and we present the analytic results as follows. If 



7* 2 > 7* + 1 (i.e., 7* > (1 + v / 5)/2), the stability condition for this strained film is 

e* 2 < 2(1 + 7*) 1 / 2 and 



*2 ^ *2 



(31) 



with ?/o 2 = (1 - ^o)/(l - 2iy o)- Otherwise, if 7* 2 < 7* + 1 (i.e., < 7* < (1 + y/E)/2), the system is stable when 



< 2 7 * 



and 



* 2 ^ # 2 
>7 > 



(32) 



2 7 * < e* 2 < 2(1 +7*) 1 / 2 and 



* 2 - * 2 



i (9 7 * 2 -2e* 4 )-2(e* 4 -3 7 * 2 ) 3 / 2 
27 7 * 2 e * 2 



(33) 



Therefore, in this stress-driven epitaxial system the instability of both morphology and composition could also occur 
above the strained spinodal temperature T° s and for large misfit e* or small solute coefficient 77*. This result is very 
different from the bulk alloy where only below the critical temperature, can the spinodal decomposition be pxesent. 
This instability is due to the coupling of morphological and compositional undulations, as pointed out by GlasE3 from 
the thermodynamic point of view. 

The corresponding stability boundary in the e*-ry* space is shown for i>q = 1/4 and different values of 7* in Fig. 
|l] (thick lines). The stability boundary for 7* = 5 corresponds to the case of 7* 2 > 7* + 1 and then is determined 
by Eq. (|3l|), while for 7* = 0.5 the conditions ([32]) and (33) are used. The stabilizing effect of surface energy can 
be seen directly from the figure, where the stable region is enlarged with the increasing value of surface tension 7*. 
The stability diagram here is symmetric with respect to the sign of misfit e* since we have assumed the composition 
independence of elastic moduli in Fig. |l|. 



B. Composition-dependent elastic moduli (E* 7^ 0, fi* 7^ 0) 

More interesting and richer results are obtained for the cases of composition-dependent elastic constants with E* 7^ 
and 7^ 0, where the coupling of misfit strain, solute strain and composition dependence of elastic moduli can highly 
affect the behaviors of perturbation growth. 

For the lattice matched films, that is, e* = with rj* arbitrary, the morphological and compositional degrees of 
freedom decouple, as obtained from dynamical equations ( pjl ) and (19). The perturbation growth rates 07, and 



also obey Eq. (j27|), corresponding to the stability properties the same as that of E* = fi* — 0. However, when lattice 
misfit exists, i.e., e* 7^ 0, the composition dependence of elastic constants leads to substantially different results. 
In the absence of atomic size difference {q* — 0), the dynamical equations for h* and </>* remain coupled, which is 
qualitatively different from the case shown in Eq. ( ]2g| ) for composition independent elastic moduli. The coupled 
perturbation growth rate a is then governed by a quadratic equation of the form similar to Eq. ( p9| ) 

<j 2 + a\a + ao = 0, 

but with different coefficient ao: 

01 = 7*fc 4 -e* 2 fc 3 + fc 2 (fc 2 ±l), 

ao = fc 2 (fc 2 ± i)( 7 *fc 4 - e * V) + [ *5 r ( U^1^ * 4fc6 - ( 34 ) 

2(1 + u )(i - v y 



When T < T° s , it is easy to show that Re(cr) > for k <C 1 and then the system is unstable, while for T > T° ff the 
film can be stable for certain values of misfit e*, as specified in the following stability conditions: If 7* 2 > 7* + 1, the 
stability occurs for 

X - x < e* 2 < 2(1 + 7*) 1 / 2 , (35) 

with 

^±Ntf . (36 ) 

X 2(l + ^)(l-^o) 2 ' [ ' 

while for 7* < 7* + 1, there are two regions of stability. In the first one, 

X- 1 < e* 2 < 2 7 *, (37) 
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is fulfilled. In the other, four conditions have to be fulfilled: 



and (iv) 



z 1 *2 ° 1 

3* < £ < 9* ' 



(i) 2 7 * <e* 2 <2(l+ 7 *) 1 / 2 , 

(ii) T * > 

(iii) 7 * > 



3(3xe* 2 /2-l)V2' 



(38) 
(39) 

(40) 



° — 1 i *2 t 

> nX , and 7 > — 



(!) 



3/2 



., I 3 *2 



e*" > gX" 1 , and 7 * 2 < — 



2 X£ " 4 



3/2 



e * - 3 ( -xe*^ - 1 



(41) 



The stability diagram of |e*| vs 7* corresponding to Eqs. (p5|)-(|4p|) is plotted in Fig. ||, where the parameters vq = 1/4, 
E* = —0.4, and p,* = —0.1 are chosen. One can see from the diagram that in the absence of compositional strain 
but with the composition dependence of elastic constants, the system above the effective critical temperature can be 
stabilized for intermediate magnitudes of misfit e* and large enough effective surface tension 7*. 

For the most general case e* 7^ 0, rf 7^ and E* 7^ 0, p,\ 7^ 0, corresponding to the lattice mismatched and 
compositionally stressed film with composition-dependent elastic constants, the coupled dynamical equations are 
described in ( |l8| ) and (|l9|), with joint perturbation growth rate a given by Eq. (^5|). The characteristic equation (pjj) 
is in fact quadratic, with coefficients 

01 = 7 *fc 4 - e* 2 fc 3 + k 2 [k 2 ± 1 + f3e*r)*] , 
a = fc 2 (7*fc 4 -e* 2 fc 3 ) [fc 2 ±l + /3eV] 
k 5 



2(l + l / )(l-^ ) 2 



(1 - 2u )e*rj* + ae* 2 2(1 - v )e*rf + ae* 2 



(42) 



with the parameters 



a = 2£*-(l + ^K, 
= 8^-5(1 + 1*)^ 



2d -o 2 ) • (43) 

The stability conditions can be derived by studying the real part of the solutions a = ^— a\ ± \/ af — 4ao^j /2, and we 

present the analytic results below for both T > T c off and T < T c cff . 
The stable epitaxial film should first fulfill: 



/3eV ± 1 > 0, 



(44) 



and then similar to the other cases, the conditions for 7 * > 7* + 1 or 7* < 7* + 1 are different. For 7* > 7* + 1, 
the stability conditions are 



e* 4 < 4(7* + l)(/?eV ± 1), 



(45) 



and 



A < 0, 



A > and 77* > 



A > and rj* < 



A 1 / 2 



4(l-H))(l-2n))' 

-pe* - A 1 / 2 
4(1-2,0X1-21/0)' 



(46) 
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where 

p = -2Ei + (l + vo)(2-voM, 

A - [p 2 - 8(1 - is )(l ~ 2is )a 2 ] e* 2 ± 16(1 + v )(l - v ) 3 (l -2u ). (47) 
On the other hand, for 7* 2 < 7* + 1 the system is stable only when the following conditions are fulfilled: 

e* 4 <4 7 * 2 (/3e*?7*±l), (48) 

as well as all the conditions in Eq. (|46|), or 

4 7 * 2 (/3e*77* ± 1) < e* 4 < 4( 7 * + l){/3e*r)* ± 1), 
and 

1 ~ 2 ^0 *2 P * * *2 .I *2 

.1-1^ + 2(1 + ^ )(1- I , ) 2e ^ 2(l + t,o)(l-^o) 2e + J" 

+^2 [97* 2 (/?e^* ± 1) - 2e* 4 ] - ^ [e* 4 - 3 7 * 2 (/? £ V ± 1)] 3/2 > 0. (49) 

Note that in Eqs. (^2|), (Q, (^H|), (^7|)-(|49|), the top sign applies when the temperature T is above the effective 
critical temperature T° s , and the bottom sign corresponds to T < r° ff . 

The stability diagrams can be calculated according to above results (^J)-(^). Here we use two sets of parameters 
(1 and 2) to plot the stability diagrams of T > Tf and T < Tf, as shown in Figs. |-|. For the first set (Set 
1, used in Figs. || and ||), where all the material parameters (e.g., T c , 7, aj, N v , and elastic constants) are chosen 
to qualitatively represent the SiGe alloy, we have vq — 1/4 and 7* = 5 obtained from Eq. (^4|), and assume that 
El = —0.4, pi = —0.1. The second set 2 is expected to qualitatively represent the InGaAs alloy and applies to Figs. 
|] and ||, where we choose vq — 1/3, 7* = 3.5, E{ = —0.25, and pi = —0.5. 

Compared to symmetric diagram Fig. for El = pi = 0, Figs. |-| exhibit the asymmetry for compressive and 
tensile layers, i.e., the stability depends on the sign of misfit e* , which is one of the major consequence of local 
composition dependence of elastic constants. In Figs. || and [|, corresponding to parameters of Set 1 (similar to 
the case of SiGe), the stabilization mainly occurs under tensile strain (e* < 0) and the stable regions increase with 
larger value of 77*, while for compressive films the instability can not be suppressed for most of the parameter values, 
especially for T < T° . In contrast, the InGaAs-like parameters of Set 2 lead to opposite asymmetry: Larger part of 
stable region is found in positive misfit e*, and layers subject to tensile strain exhibit less stability, as shown in Figs. 
I and § 

The other important effect of composition-dependent elastic moduli is on the system below effective critical temper- 
ature . For all the other cases described above, including the ones for El ^ and pi ^ but with one of e* and 77* 
equal to zero, the compositional profiles for T < Tj? are unstable, in agreement with the usual expectation that the 
strained alloy near 50 — 50 mixture should exhibit decomposition and phase segregation below the coherent spinodal 
temperature T° s . However, the coupling of all the factors of misfit strain, solute strain and composition-dependent 
moduli causes different and new effects. As shown in Figs. || and ^ when all the variables e*, 77*, El and pi are 
nonzero, the film below Tf can also be stable for certain range of parameters. That is, it is possible to suppress the 
surface decomposition even for T < T° s due to the coupling effects in this heteroepitaxial system. 



C. Oscillatory instability (Im(<r) 7^ 0) 

When the system corresponds to the unstable parameters region of stability diagram, that is, Re(er) > 0, the 
imaginary part of a determines whether the onset of this instability is steady (Im(cr) = 0) or oscillatory (Im(cr) ^ OX 
The occurrence of oscillatory instability has. been found in the study of directional solidification for stressed solidlEl 
and the growing process of alloy thin films ,E3 and been attributed to the phase difference between surface morphology 
and composition field, induced by nonlocal elastic stresses. 

We calculate the imaginary part of perturbation growth rate a through characteristic equation ( |25| ) . If for a certain 
wavenumber k, we have both Re(cr) > and Im(cr) 7^ 0, the oscillatory instability may occur. The results for e* ^ 0, 
77* 0, and T > T° s , corresponding to the parameters range of most experiments on strained films, are shown in Figs, 
[j], and |i| with the thin dashed or dotted boundary curves. For the case of composition independent elastic moduli, 
as shown in Fig. [I], the oscillatory unstable regions are symmetric with respect to misfit and more regular. Large 
solute coefficient 77* favours the occurrence of oscillatory instability, and if misfit e* is large enough, the oscillatory 
instability is obtained when 77* exceeds a fixed value ?7g = [(1 — v^)/{\ — 2^o)] 1 ^ 2 , which can be derived analytically. 
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When the elastic constants are composition dependent, i.e. E* ^ and /j* ^ 0, the oscillatory regions are asymmetric 
and irregular (see Figs. || and ||). Note that in these parameter regions of oscillatory instability, steady instability can 
also exist (but for different mode k), and the competition between these two kinds of unstable modes determines the 
surface profile. When oscillatory modes dominate or coexist with steady modes, the surface disturbance is expected, 
to propagate laterally, with the phenomenon that one side of the surface bump will grow faster than the other side.tS 

V. DISCUSSION AND CONCLUSIONS 

Our calculations above have shown that the stability problem of free alloy film surface under both morphological 
and compositional strains is essentially a nonequilibrium and dynamical problem even for static films, i.e., in the 
absence of growth. Although in the sense of thermodynamics and equilibrium, it has been demonstrated that the 
instability can appear in any stressed alloy with free surface,E3 the physically based choice of non-equilibrium evolution 
dynamics leads to different conclusion that the system can be stabilized for certain values of parameters. As shown 
in Figs. the joint instability can be suppressed by large enough compositional strain (r)*), which has also been 
found for growing filmsjaa the film- vapor local equilibrium model of Guyer and VoorheesE3 as well as in the dynamical 
model of Spencer et al. Hjsrith unequal atomic mobilities for different alloy constituents. 

In the previous worloOo this stabilization of compositional stresses would be overcome by larger magnitude of 
misfit e*, while in our results similar phenomenon occurs for e* values far from zero, but the other parts of stability 
diagrams are more complicated. In Fig. |] (with parameters of Set 2 and T > T° ff ) and Figs. || and |] for T < T° s , 
when the magnitude of misfit e* is close to the minimum of the stability boundary and is made smaller, then higher 
value of rj* (related to larger compositional strain) is needed to suppress the instability, and for parameters of Set 1 
with T > Tf (Fig. g) the stable regions are much more irregular, due to the combination of misfit and compositional 
strains as well as composition-dependent elastic moduli. Even for the case of composition-independent moduli, the 
results herje-|^hown in Fig. |l|) are different from before. For sufficiently large 77, the nonphy-sical short wavelength 
divergenceE2lE3 is avoided due to the inclusion of gradient energy and the return of instabilitycll is not found here due 
to the consideration of coupling between surface and bulk phases. 

The introduction of composition dependence of alloy elastic constants, which makes the effective elastic effects 
nonlocal, leads to the presence of asymmetry in the stability of films under compression or tension. Some experiments 
have tested this misfit sign dependence. Although these observations are all for the growing films, they can still be 
helpful to check our theoretical results for films without deposition. Our calculations using parameters similar to that 
of SiGe (shown in Figs. || and ||) exhibit t he p pr a fa rence of stability for tensile layers, consistent with the experimental 
findings of Xie et al.u Various experiments" 3 H 14 ti 15 l indicated that the way of asymmetry for InGaAs alloys depends on 
the materials deposition rate, while our results in Figs. [| and || with a specific selection of parameters (Set 2) suggest 
that the compressive layers are more stable for non-growing film. 

The simultaneous interaction of misfit strain, solute strain and composition-dependent clastic constants makes the 
stability possible even for T < T° ff , which is not possible in the absence of any one of them. The corresponding 
diagrams are shown in Figs. || and |^. This phenomenon is very different from what we expect, since in usual bulk 
strained alloy the spinodal decomposition always occurs below T^ S S3 The related experiments are lacking since most 
of the epitaxial experiments are carried out above the effective critical temperature T° ff . 

Note that although in the above analysis, we have distinguished the stability results and diagrams into two tem- 
perature regimes T > T° s and T < T^ s , in each regime the stability properties are still temperature dependent. This 
can be seen from our theoretical diagrams (Figs. |l| - ^|) as well as the charateristic equation (|2^) for perturbation 
growth rate, where the system stability is shown to depend on the rescalcd parameters e*, r/* , and 7*, which in fact 
are all proportional to | 1 1 / 2 with r linearly dependent on temperature T (see Eqs. (|24|) and (p3[)). This temperature 
dependence of stability property can also be obtained from Eq. ( p6| ) that has dimension-full form. 

It is interesting to compare our theoretical results for instability with the.-Qbse£vations of isothermal annealing 
experiments on strained films. Experiments on SiGe/Si postdeposition systeml3'E9 : Eil have exhibited a morphological 
evolution procedure from an initial planar film surface to a rough surface profile with ripples or islands during 
annealing. Also, along with this morphological modulation, the Ge segregation has been found on the surface of 
SiGe layers by Walther et al.p which corresponds to the coupling of compositional and morphological instabilities 
studied here. The temperature effect on the stability of Sio.sGeo.s strained films has been investigated by Chen et 
aZ.Ej They observed that the evolution to unstable surface morphology only occurs above an annealing temperature, 
and attributed this sharp temperature dependence to the existence of an energy barrier. Our theoretical results (for 
the case of T > T£ s , e.g., stability diagram like Fig. ||), obtained from the surface diffusion mechanism, can explain 
this temperature dependent phenomenon without the introduction of an energy barrier. For small temperature T 
(but still above T° s , as in real experiments), we have small \r\ (see Eq. (p3[)) and then large 7* (resulting in the large 
stable region of the stability diagram) and large 77*, which are apt to suppress the instability. When T increases, 
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the value of e* becomes smaller, and more importantly, 7* and rf also decrease, rendering the system closer to the 
unstable region of the stability diagram. Thus, one can reach a transition temperature above which the annealing 
system is unstable, as found in the experiments. 

In summary, we have developed a continuum model to study the nonequilibrium evolution processes of strained 
alloy films in the absence of growth. With the consideration of coupling between surface and bulk states as well as the 
composition dependence of elastic moduli (both E and /i), new and more complicated stability results and diagrams 
have been obtained using linear stability analysis. In general case, joint stability or instability is found due to the 
coupling between morphological and compositional perturbations. More importantly, the interplay of morphological 
and compositional strains as well as the composition-dependent elastic constants leads to the stability dependence 
on the sign of film-substrate misfit, and the possibility of stabilization even for films below T° s . Here we only study 
the early film evolution, and the nonlinear effects should be considered for later regimes and for the determination of 
detailed surface morphologies and patterns. 
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FIG. 1: Stability diagram for non-growing strained alloy film with E'l — fi\ = 0, temperature T > T° s , and vq = 1/4. 
Region marked by "S" is the stable region. Long dashed and solid thick lines denote the stability boundaries for 7* = 0.5 
and 5, respectively. In the unstable region, the domain marked as "Ug" corresponds to steady instability and that marked as 
"^S&o" corres P on ds to steady or oscillatory instabilities depending on the wavenumber k. The boundaries between "Ug" and 
"^S&O" re gi° ns are indicated by dashed (7* = 0.5) and dotted (7* = 5) thin curves. 



T>T C , r|*=0, d =1/4, E^-0.4, m*=-0.1 



U 




1.5 1.6 1.7 



1.8 1.9 2.0 
|e| 



2.1 2.2 2.3 



FIG. 2: Stability diagram for 77* = non-growing strained alloy film, with T > T c cff , v = 1/4, E{ = -0.4, and n\ = -0.1. 
Stable and unstable regions are marked as "S" and "U" , respectively. 




FIG. 3: Stability diagram for non-growing strained alloy film with T > Tj? . Parameters of Set 1 are chosen: — —0.4, 
fil — —0.1, uo — 1/4, and 7* = 5. Stable and unstable regions are marked in a manner similar to that in Fig. pi 




FIG. 4: Stability diagram for T > T° strained alloy film in the absence of growth, with Set 2 parameters: E{ — —0.25, 
Hi = —0.5, uo — 1/3, and 7* = 3.5. Regions with different stability properties are marked as in Figs. □ and H. 
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FIG. 5: Stability diagram for non-growing strained alloy film with T < T° and composition-dependent elastic moduli. The 
parameters are of Set 1 as described in Fig. H. 
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FIG. 6: Stability diagram for non-growing strained alloy film with T < T£ as well as the Set 2 parameters as described in Fig. 
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